Dark Matter in the Solar System 
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We determine the density and mass distribution of dark matter within our Solar System. We 
explore the three-body interactions between dark matter particles, the Sun, and the planets to 
compute the amount of dark matter gravitationally captured over the lifetime of the Solar System. 
We provide an analytical framework for performing these calculations and detail our numerical 
simulations accordingly. We find that the local density of dark matter is enhanced by between three 
and five orders of magnitude over the background halo density, dependent on the radial distance 
from the Sun. This has profound implications for terrestrial direct dark matter detection searches. 
We also discuss our results in the context of gravitational signatures, including existing constraints, 
and find that dark matter captured in this fashion is not responsible for the Pioneer anomaly. We 
conclude that dark matter appears to, overall, play a much more important role in our Solar System 
than previously thought. 
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I. INTRODUCTION 

Accounting for about 23% of the energy density of the 
Universe dark matter is an integral part of our sur- 
roundings. It dominates the more familiar, baryonic com- 
ponents, which comprises only 4.4% of the Universe, on 
the largest scales. Achieving an understanding of this 
perplexing dark element is one of the most compelling 
unsolved problems in modern astrophysics. 

The astrophysical evidence for the existence of dark 
matter is overwhelming ,as observations of the cosmic 
microwave background the power spectrum of the 
Universe [|J , and colliding galaxy clusters [l| all point to- 
wards the same conclusion. With an understanding that 
a dark, pressureless, fluid-like component permeates the 
Universe, astrophysicists have successfully simulated the 
large-scale processes of structure formation in the context 
of a ACDM Universe Q. More recently, some attention 
has been given to the dark matter substructure formed 



on subgalactic scales, down to scales of order 10 pc 
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However, relatively little consideration has been given to 
the distribution of dark matter within our own Solar Sys- 
tem. 

Yet, dark matter may prove to be profoundly impor- 
tant in our Solar System for both its additional gravita- 
tional effects on planets and other orbiting bodies [1,0, 
as well as the motions of spacecraft (IHllO). Furthermore, 
a knowledge of the density and velocity of dark matter 
particles is particularl y im portant for terrestrial direct 
detection experiments fllj . 

In this paper, we model the Solar System and the dark 
matter that it encounters in order to quantify how much 
dark matter we expect the Solar System to have captured 
over its lifetime. Through favorable three-body gravita- 
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tional interactions between a dark matter particle, the 
Sun, and any of the planets, a non-zero and possibly 
significant fraction of the dark matter passing through 
the Solar System will become gravitationally bound to 
it. The remainder of this paper is focused on solving 
this problem, and is laid out as follows: section 2 details 
the model chosen for the Solar System and galactic prop- 
erties, and provides analytic details of our calculations. 
Section 3 sets forth the computations undertaken to suc- 
cessfully determine the probability of binding dark mat- 
ter to the Solar System. Section 4 presents our results, 
including the density and mass distributions of dark mat- 
ter with respect to distance from the Sun. We recom- 
mend that anyone not interested in the details of our 
calculations skip directly to section 4. Finally, section 
5 concludes with a discussion of the conclusions reached 
from our analysis, detailing significant implications for 
dark matter detection and presenting a comparison of 
our results with the current experimental and observa- 
tional constraints. 



II. DARK MATTER CAPTURE 

Our Sun (and hence our Solar System) is presently 
moving through the galaxy with well-known parameters 
[l2T | that have changed little, despite refinements in mea- 
surements [13j . over many years. More recently, we have 
been able to determine that our Milky Way, like all com- 
parable galaxies, is pervaded by a dark matter halo with a 
specific density profile [3] . A-body simulations also pro- 
vide insights into modeling the Milky Way [l5j]. Based 
on the fact that our Sun is not an isolated body, but 
rather has eight planets as well as a number of other, less 
significant bodies gravitationally bound to it, a non-zero 
fraction of this dark matter will be captured by favorable 
three-body interactions between a dark matter particle, 
the Sun, and an orbiting body. 

In order to compute the amount of dark matter cap- 
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tured by the Solar System over the course of its 4.5 Gyr 
history we assume the galactic parameters shown in Ta- 
ble U below. We then find that the Sun moves through 



Parameter 


Value 


Oort's Constant A (kms~ kpc - ) 


14.4 


Solar Radius r Q (kpc) 


8.0 


Speed of Local Standard of Rest (LSR) v\ sr 


200 


Sun's speed relative to LSR W0,i sr 


20 


Velocity relative to the Galactic Plane v z 


7 


Solar Period relative to the Galactic Plane t osc (Myr) 


63 


Local Dark Matter density Pdm(t"©) (Mq pc -3 ) 


0.009 


Mass within tq of Galactic center M cnc (rQ) (10 10 Mq) 


9 


Local rms— velocity of Dark Matter u rms (r0) 


220 



circularized velocity around the Sun v p . (This analysis 
will be repeated eight times, once for each of the plan- 
ets.) We first choose a very large r M to be the vector 
distance from the dark matter particle to the Sun, so 
that | Too | 3> \f p \, but small enough that when we com- 
pute the dark matter particle's angular momentum with 
respect to the Sun, we get a reasonable (i.e., non- infinite) 
value. 

We then perform a coordinate transformation to shift 
into the Sun's rest frame, computing the velocity of the 
dark matter at infinity (i>oo) in that frame: 



("DM - Vq) ■ r a 



("DM —Vq) X r c 



(3) 



TABLE I: The local galactic and dark matter parameters used 
for the Milky Way and the present values of the Sun's distance 
and velocity components. Unless otherwise noted, all veloci- 
ties are in units of kms -1 . 



the galaxy with a velocity (i>©) with components given 
in cylindrical coordinates: 



v GtP = Ar Q cos (21), 
vq,4, = Uisr - u ,isr sin (21), 
'2ttA 



v z sin | - — I , 

( £osc / 



(1) 



where / is the galactic longitude, t is a time coordinate, 
and the remaining parameters are defined in Table|TJ The 
individual dark matter particles are assumed to follow a 
Maxwell-Boltzmann distribution fl6| with a probability 
distribution function f(v) given by 



f(v) 




(2) 



where the local rms- velocity w r ms('*©) is given in Table 
III 

Therefore, the Sun sweeps out a predictable three- 
dimensional path over its 4.5 Gyr history. With the dark 
matter having a local density Pdm(Tq), a rms-velocity 
VrmsC^o) and a velocity distribution as given above, the 
fraction of dark matter captured can be calculated in a 
straightforward fashion via modeling of the Solar System 
and the dark matter particles passing through it. How- 
ever, the number of dark matter particles encountered is 
far too large and the rate of particle capture is far too 
small to effectively simulate using TV-body methods. We 
are therefore forced to use analytic approximations to 
shape this problem into a tractable one. 

We begin by considering an ensemble of dark matter 
particles at infinity each with a speed |udm| chosen from 
the Maxwell-Boltzmann distribution and a random ori- 
entation vum- We also consider the Sun, and a single 
planet with mass m„, distance from the Sun r p , and a 



where v^n and Woo_l are the components of the dark mat- 
ter particle's velocity parallel and perpendicular to the 
Sun's, respectively. 

The dark matter particles we are interested in, as far 
as the possibility of gravitational capture goes, are those 
that pass within a distance r p of the Sun. The dark mat- 
ter particles, at infinity, will have an angular momentum 
L given by 



mv 00 ±r 



(4) 



where m is the mass of a dark matter particle. We note 
that this mass is completely unimportant in our analy- 
ses, as it does not enter into any of our equations; only 
the combination L/m appears. These particles are all 
in hyperbolic orbits around the Sun initially, with eccen- 
tricities e given by 




G 2 Ml 



(5) 



We then find that any particle that meets the following 
condition will pass within a distance r p of the Sun: 



e > 



(L/m) 



GM, 



1 



(6) 



Upon encountering the planet, the particle may receive 
boosts (or anti-boosts) at two points during the interac- 
tion. The first occurs at entry into the sphere of radius 
r p and the second occurs upon exit of the sphere. Taking 
both of these opportunities into account, we consider the 
approximation that the planet's position at any time is 
given by a random location on the sphere of radius r p . 
We determine that of the particles that pass through the 
sphere of radius r p , a fraction r^/r 2 of those will gravi- 
tationally encounter the planet, where equals 



n = 1.1 5r p 



Mr. 



1/3 



(7) 



with Tb defined to be the radius of a "sphere of influence" 
of a smaller gravitational body relative to a larger one 
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FIG. 1: Geometric setup for a dark matter particle approach- 
ing the Solar System from infinity in the Sun's rest frame. The 
semi-major axis of the planet in question's orbit is shown by 
r p , with the velocity that the particle strikes the imaginary 
sphere of radius r v given by v ln - Through the conservation of 
angular momentum, the angle 9 can be determined, as shown 
in equation © , and the angle that the velocity v ln makes with 
r p , defined as (f>, is given by equation (|10[) . 



We then need to determine whether the gravitational 
encounter is favorable enough to transition the dark mat- 
ter particle from an unbound, hyperbolic orbit around 
the Sun to a bound, elliptical one. We begin by deter- 
mining the velocity of the incoming dark matter particle 
at a distance r p from the Sun (= Ui n ). The geometry of 
an incoming dark matter particle as it enters the Solar 
System, possibly encountering a given planet's sphere of 
influence, is illustrated in Figure [1] This velocity will 
have a magnitude given by the conservation of energy, 



IGMr. 



(8) 



and we are interested in determining the components par- 
allel to the radial vector towards the Sun (i>i n ||) and per- 
pendicular to it (v- m ±). We can determine these, in the 
Sun's rest frame, by first finding the angle of deflection 
(6) caused by the gravitational force on the position of 
the dark matter particle's azimuthal angle from infinity 
to r„, 



(L/r, 



GMqTj, 



1 



(9) 



and then, through the conservation of angular momen- 
tum, the angle of deflection (0) of the particle's velocity, 



L/m 



sin(0) ±cos(6») 1 



{L/ m y 



( 10 ) 

where 4> is chosen to be the smaller in magnitude of the 
two possible angles. From equations ([5] fTU|) . we can then 



determine the components of ^i n to be 

w in|| = Wm [sin(0) cos(0) — cos(#) sin(0)] |, (11) 



I'i, 



m|| 



(12) 



Assuming that the particle does get within a distance 
r b of the planet (which, as stated above, it does with 
probability r^/r^), we then need to determine the gravi- 
tational effect of the planet on the dark matter particle's 
orbit with respect to the Sun. We assume that the planet 
is positioned randomly in space at a distance r p from the 
Sun, and moves with a speed \v p \ in a random direction 
perpendicular to v- ln \\. 

We then switch to the planet's rest frame, and obtain, 
for the dark matter particle, a velocity (t>dm) with com- 
ponents 



«in||) 

«in_L - WpSin(a), 
—Vp cos(a), 



(13) 
(14) 
(15) 



where a is a random angle between and 2tt, and a po- 
sition with respect to the planet (rdm) with components 



Tdm — {\Jrl - r\ m y - rj m z ,r dm ^,r dmjZ }, (16) 

where fd m ,y and rd m ,z are randomly chosen to lie within 
a circle of radius r b in the yz— plane. We also note that in 
this rest frame, the Sun appears to move with a velocity 



vq = (0, — UpSin(a), — w p cos(a)). 



(17) 



We then compute the components of i>d m parallel and 
perpendicular to rd m , 



Wdm|| 



Vdm± 



Vdra X Tdm 



n 



(18) 
(19) 



and the components of vq with respect to these new di- 
rections, 



v &\\ 

VQ3d 



Vq ■ r dm 



■ 



"dm X Tdm 
|Udm_L r b \ 



V Q V Q\\ U 0-L' 



(20) 
(21) 
(22) 



where «03d is the component of w© orthogonal to both 
Vdm\\ and Wdmi- In this coordinate system, it is easy to 
compute the angle (/3) that the planet causes the dark 
matter particle to deflect by, 



(3 = 2tan~ 1 



Gm 



n Vdm± "dm 



(23) 
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where the final outgoing velocity of the dark matter par- 
ticle (vout) is then given by 



Vout = (vdm\\ cos(P) + v dm j_ sin(/3), 

-Udm|| sin(/3) + w dm j_ cos(/3), 0). 



(24) 



At last, we can compute the total final speed of the dark 
matter particle (|v/|), as it leaves the sphere of influence 
of the planet, in the rest frame of the Sun, 



\Vf\ = \J (Vout - v Q ) ■ (u out - Vq). (25) 

The particle will be gravitationally captured if 



Vf <C ^escape ypj 



'2GMr. 



(26) 



and become bound in an elliptical orbit about the Sun 
with semi-major axis 



= / \v f \ 2 2 
\GM Q r p 



(27) 



A point to note is that we do not consider any further 
three-body interactions between captured dark matter, 
the planets and the Sun. Once the dark matter is cap- 
tured, we assume it remains captured without any further 
gravitational interactions of significance. However, such 
a back-reaction will exist, and could potentially decrease 
the amount of dark matter bound to the Solar System 
by a significant amount. 



III. COMPUTATIONS 

For each of the eight planets in our Solar System, we 
perform the calculations outlined in section[n] The plan- 
ets are assumed to have the parameters shown below 
in Table [Til We begin by creating a cumulative distri- 



Planet 


Distance to Sun 


Mass 


Speed 


Mercury 


5.79 x 10 7 km 


1.68 x 10~ 7 M Q 


48 kms -1 


Venus 


1.08 x 10 s km 


2.46 x 10~ 6 M Q 


35 kms - 


Earth 


1.496 x 10 8 km 


2.99 x 10~ 6 M© 


30 kms -1 


Mars 


2.28 x 10 s km 


3.21 x 10~ 7 M Q 


24 kms -1 


Jupiter 


7.78 x 10 8 km 


9.50 x 10" 4 Mq 


13kms _1 


Saturn 


1.43 x 10 9 km 


2.86 x 10" 4 M Q 


9.6kms _1 


Uranus 


2.87 x 10 9 km 


4.40 x 10" 5 M Q 


6.8 km s -1 


Neptune 


4.50 x 10 9 km 


5.11 x 10~ 5 Mq 


5.4kms _1 



TABLE II: Planet-Sun distances, planetary masses and 
speeds for the eight Solar System planets considered in this 
analysis. Distances are given in units of km, masses are given 
in units of solar masses, where Mq = 1.9884 x 10 30 kg, and 
speeds in units of kms -1 . 



bution function for the Maxwell-Boltzmann distribution 



(equation ([2])), obtaining cumulative probability P(v) as 
a function of the dark matter's velocity, 



P(v) = f f(v')dv' 
Jo 




where P(v) is the probability of finding a dark matter 
particle with velocity less than or equal to v. From this 
cumulative distribution function, we determine the ve- 
locity with respect to the Sun and keep only those dark 
matter particles that meet the condition 



< 2 



\2GMr. 



1/2 



(29) 



where v p and r p are the values for the appropriate planet 
as given in Table HT1 and Voo is the dark matter particle's 
speed with respect to the Sun's reference frame at in- 
finity. We choose the condition in equation ([29]) because 
even the most favorable gravitational interaction with the 
planet can only decrease the speed of an incoming dark 
matter particle by 2v p . The condition in equation (|29|) 
is such that a particle reaching a distance r p from the 
Sun will have a velocity no greater than 2v p . We gener- 
ate at least one million unique particles that satisfy this 
condition for each planet. 

Keeping track of the total number of particles simu- 
lated before generating the one million we seek, we com- 
pute the probability of particles satisfying the preceding 
speed constraint. We then generate random directions 
for a particle, with an initial position at r^, and a speed 
chosen randomly from the the set of one million parti- 
cles. We simulate a total of 100 billion particles for each 
of the Jovian planets and from 300 to 600 billion particles 
for each of the inner, rocky planets, which require more 
due to their much smaller distances from the Sun. We 
adopt ^ = 2.67 x 10 15 m, placing it orders of magni- 
tude beyond the orbit of Neptune, but still close enough 
so that a reasonable number of simulated particles will 
interact with each planet, based on the condition in equa- 
tion ([SI). We then compute the location and velocity of 
all of the particles that do pass within a distance r p of 
the Sun (correcting for the fact that the Too chosen is not 
actually infinite), as well as the probability of randomly 
generated particles passing within r p of the Sun. 

We then reduce that probability further by a factor of 
rl/r p , as discussed in section [ill as only that fraction of 
the dark matter particles passing within a distance r p of 
the Sun will pass within a distance r\, of the planet in 
question. 

We then generate, from the particles that have passed 
all the cuts up until now, a random set of positions for the 
planet within a three-dimensional distance r\> of the dark 
matter particle, according to equation ([1~6]) . along with 
velocity directions for the planet in accordance with equa- 
tions {131-15). From the particles that survive the earlier 
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cuts, we choose enough random positions and velocity 
directions to generate one billion particles for this final 
step. By boosting to the planet's rest frame, calculating 
the change in direction due to the gravitational interac- 
tion between the dark matter particle and the planet, and 
then boosting back to the Sun's rest frame, we obtain the 
final speed of the particle. We tabulate the particles that 
become bound to the Sun as a result of this final inter- 
action, and compute both the total probability of grav- 
itational capture and the distribution of the semi-major 
axes of the captured particles. 



IV. RESULTS 

Given the chosen in section IIIII and the fact that 
the Solar System has had approximately 4.5 billion years 
(defined as the lifetime of the Solar System, igg) to accrue 
dark matter via this mechanism, we determine that the 
total amount of dark matter encountered is given by 



M DM = pvM{r & )V 

= PDM(r©)7rr^«0tss - 203 M Q , 



(30) 



where Vq is the average velocity of the Sun with respect 
to the galaxy, determined to be 208 km s -1 using a time 
average of the data from equation |T]) and Table fl] Our 



Planet 


Fraction Captured 


Mass Captured 


Mercury 


1.03 x KT 16 


2.09 x 10~ 14 M© 


Venus 


8.71 x KT 16 


1.77 x 10~ 13 M© 


Earth 


9.41 x 10~ 16 


1.91 x 10~ 13 M© 


Mars 


2.91 x lCT 16 


5.91 x 10 -14 M© 


Jupiter 


1.23 x 10~ 13 


2.50 x 10~ 14 M© 


Saturn 


7.06 x 10~ 14 


1.43 x 10~ n M© 


Uranus 


2.87 x 10" 14 


5.83 x 10~ 12 M© 


Neptune 


3.98 x 10~ 14 


8.08 x 10" 12 M© 



TABLE III: Fraction and total mass of all dark matter par- 
ticles captured due to gravitational interactions with each of 
the eight planets. Although the absolute numbers are much 
smaller for the inner, rocky planets compared to the Jovians, 
they are still significant for determining the densities of dark 
matter, as they dominate at radii smaller than half of Jupiter's 
semi-major axis. 



calculations and computations then allow us to determine 
what fraction of this total mass winds up getting gravita- 
tionally captured via these three-body interactions. Our 
results are presented in Table IIIII 

The two most important factors contributing to cap- 
ture are the mass and the orbital semi-major axis of the 
planet. The more massive and the further away a planet 
is from the Sun, the more effective it will be at capturing 
dark matter, as the cross-section for favorable gravita- 
tional interactions is given by the radius of the sphere of 
influence (r&) from equation ([7]). 
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Top: DM captured by Solar System 
Bottom: Background Galactic DM 
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FIG. 2: Dark matter density vs. distance from the Sun in our 
Solar System at the present day. Density is presented in units 
of kg AU~ 3 , with distances given in units of AU. 



Equation (j2"T|) allows us to obtain a measure of the spa- 
tial distribution of the captured dark matter particles. In 
Figure [H we illustrate the local density of dark matter, 
taking the semi-major axis of the captured dark matter 
particle as a proxy for its position, as a function of dis- 
tance from the Sun. For comparison, the background 
galactic halo density is also shown. Of particular note is 
the dark matter density at 1 AU, which is greater than 
the background dark matter density (from the underlying 
galactic halo) by a factor of 1.63 x 10 4 . 

Figure [3] indicates the total mass of dark matter en- 
closed within a certain radius from the Sun due to both 
the captured dark matter and the underlying galactic 
halo. Between approximately 0.2 AU and 100 AU, the 
amount of dark matter bound to the Solar System is 
much more massive and dense than the background halo. 
Within the orbit of Mercury, Earth, Mars, and Neptune, 
we find that approximately 1.91 x 10 16 kg, 3.23 x 10 17 kg, 
4.87 x 10 17 kg, and 7.69 x 10 19 kg of dark matter is en- 
closed, respectively. The total amount of matter dark 
bound to the Solar System as the result of gravitational 
capture is 1.07 x 10 20 kg, or 1.78 x 10~ 5 Earth masses. 



V. DISCUSSION 

We predict the presence of a new component of dark 
matter within the Solar System due to gravitational 
capture. We find that, within the orbit of Neptune, 
7.69 x 10 19 kg of dark matter has become bound to our 
Solar System due to the capture mechanism over its life- 
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Distance from Sun (AU) 



FIG. 3: Cumulatively enclosed mass as a function of distance 
from the Sun at present. This is representative of the addi- 
tional mass gravitationally felt by an object in orbit around 
the Sun due to the presence of captured dark matter. 



time. This is about a factor of 300 greater than the 
background mass from the galactic halo, as shown in 
Figure O The density of the captured dark matter is 
also significantly enhanced in comparison with the halo, 
as demonstrated in Figure [2l At the Earth's orbital ra- 
dius, density is enhanced by more than four orders of 
magnitude over the local halo density, with a value of 
3.3 x 10 16 kgAU~ 3 . Such elevated levels of dark matter 
have important implications for direct detection experi- 
ments and can be tested as a potential explanation for 
spacecraft anomalies. 

Direct detection searches for dark matter rely predom- 
inantly on nuclear recoils [111 ], where the rate of inter- 
action is dependent on the dark matter's density, veloc- 
ity, and interaction cross-section (which itself may have 
a velocity dependence). Our determination of the lo- 
cal dark matter density and velocity distribution are of 
great importance to direct detection experiments. The 
most recent calculations that have been carried out as- 
sume that the properties of dark matter at the Sun's lo- 
cation are derived directly from the galactic halo [H[ . By 
comparison, we find that terrestrial experiments should 
also consider a component of dark matter with a density 
16,000 times greater than the background halo density, 
albeit with a v lms on the order of the Earth's orbital 
speed (30kms _1 ), about a factor of seven smaller than 
the rms-velocity of the unbound halo particles. If this 
dark matter is not efficiently ejected by future interac- 
tions after the initial capture, the gravitationally bound 
component of dark matter may wind up dominating the 



signal expected from future direct detection experiments. 

One method that has been used to constrain the 
amount of dark matter in the Solar System has been 
careful, long-term observations of planetary and satel- 
lite motions. Constraints have been obtained both from 
planetary orbital data [6j, |8[ and perihelion shift obser- 
vations [3, [l9| . The most stringent results for the dark 
matter density near Earth constrain £dm to be less than 
6.0 x 10 16 kg AFT 3 from orbital data [|] . For the densities 
near Mercury and Mars, perihelion precession provides 
the tightest constraints, yielding upper limits on /?dm of 
8.7 x 10 18 kg AU" 3 near Mercury and 5.4 x 10 16 kg AU" 3 
near Mars fl9j| . Our results satisfy these constraints, 
as we predict the densities near Mercury, Earth, and 
Mars to be 1.5 x 10 17 kgAU" 3 , 3.3 x 10 16 kgAU" 3 , and 
8.5 x 10 15 kgAU~ 3 , respectively. Predictions about the 
effects of dark matter on planetary orbits are potentially 
observable, as the constraints on Earth and our predic- 
tions differ by less than a factor of two. 

Another interesting issue to address is the Pioneer 
anomaly and the possibility that it has arisen due to the 
dark matter bound to our Solar System. Measurements 
show that Pioneer 10 and 11 have exhibited extra accel- 
erations towards the Sun of 8.09 ±0.20 x 10~ 8 cms~ 2 and 
8.56 ±0.15 x 10" 8 cms" 2 , respectively In order for 
dark matter to have caused this, at least 3 x 10 _4 Mq of 
dark matter is required within 50 AU of the Sun Q . The 
lower bound on the dark matter density capable of caus- 
ing the anomalies is 6.0 x 10 18 kgAU -3 for an inelastic 
scattering of dark matter particles [lfj. Our results do 
not match either of these predictions, as we predict only 
~ 10 20 kg of dark matter enclosed within the entire So- 
lar System and a dark matter density that never exceeds 
2.0 x 10 17 kgAU~ 3 anywhere. We conclude that the Pio- 
neer anamoly cannot be caused by dark matter that has 
been captured by our Solar System. 

Overall, we find that dark matter in our Solar Sys- 
tem is far more important than previously thought. Due 
to gravitational three-body interactions between dark 
matter particles, the Sun, and the planets, a significant 
amount of dark matter winds up gravitationally bound 
to our Solar System, resulting in density enhancements 
between two and five orders of magnitude, depending on 
the distance from the Sun. A future direction for this 
work is to include back-reaction effects, such as subse- 
quent gravitational scatterings of the captured particles. 
These may prove to be important in ejecting a portion of 
the captured dark matter particles, and in reducing the 
total amount of dark matter that remains bound to our 
Solar System. More accurate modeling of our galaxy's 
dark matter halo, such as the possible inclusion of a dark 
matter disk in the plane of our galaxy [20] , will also alter 
the net amount of dark matter captured, and is worth 
further study. Our results may lead to exciting new di- 
rections in direct detection experiments and our under- 
standing of dark matter on the smallest scales. 



7 



Acknowledgments edges NASA grants NNX07AH11G and NNX07AC51G 

as well as NSF AST-0707725 for support. 

We thank Romeel Dave and Daniel Eisenstein for use- 
ful conversations concerning this work. E.R.S. acknowl- 



[1] E. Komatsu et al. [WMAP Collaboration], 
larXiv:0803.0547l [astro-ph], 

[2] S. Cole et al. [The 2dFGRS Collaboration], Mon. Not. 
Roy. Astron. Soc. 362, 505 (2005); W. J. Percival et al, 
Astrophys. J. 657, 645 (2007). 

[3] D. Clowe, M. Bradac, A. H. Gonzalez, M. Markevitch, 
S. W. Randall, C. Jones and D. Zaritsky, Astrophys. J. 
648, L109 (2006). 

[4] J. S. Bullock et al, Mon. Not. Roy. Astron. Soc. 321, 
559 (2001). B. Allgood, R. A. Flores, J. R. Primack, 
A. V. Kravtsov, R. H. Wechsler, A. Faltenbacher and 
J. S. Bullock, Mon. Not. Roy. Astron. Soc. 367, 1781 
(2006); P. Bett, V. Eke, C. S. Frenk, A. Jenkins, J. Helly 
and J. Navarro, Mon. Not. Roy. Astron. Soc. 376, 
215 (2007); O. Hahn, C. Porciani, C. M. Carollo and 
A. Dekel, Mon. Not. Roy. Astron. Soc. 375, 489 (2007); 

A. V. Maccio', A. A. Dutton, F. C. van den Bosch, 

B. Moore, D. Potter and J. Stadel, Mon. Not. Roy. As- 
tron. Soc. 378, 55 (2007); J. D. Cohn and M. J. White, 
Mon. Not. Roy. Astron. Soc. 385, 2025 (2008); 

[5] A. M. Green, S. Hofmann and D. J. Schwarz, JCAP 
0508, 003 (2005); A. Loeb and M. Zaldarriaga, 
Phys. Rev. D 71, 103520 (2005); E. Bertschinger, 
Ph ys. Rev. D 74, 063509 (2006); E. R. Siegel et 
al, |arXiv:astro-ph/0611864] T. Goerdt, O. Y. Gnedin, 
B. Moore, J. Diemand and J. Stadel, Mon. Not. Roy. As- 
tron. Soc. 375, 191 (2007); K. M. Zurek, C. J. Hogan and 
T. R. Quinn, P hys. Rev. D 75, 043511 (2007); L. D. Duffy 
and P. Sikivie. larXiv:0805.4556l [astro-ph]. 

[6] J. D. Anderson, E. L. Lau, A. H. Taylor, D. A. Di- 
ms, D. C. Teplitz and V. L. Teplitz, Astrophys. J. 342, 
539 (1989); J. D. Anderson, E. L. Lau, T. P. Krisher, 
D. A. Dicus, D. C. Rosenbaum and V. L. Teplitz, As- 
trophys. J. 448, 885 (1995); L. Iorio, JCAP 0605, 002 
(2006); M. Sereno and Ph. Jetzer, Mon. Not. Roy. As- 
tron. Soc. 371, 626 (2006). 



[7] O. Gron and H. H. Soleng, Astrophys. J. 456, 445 (1996); 
I. B. Khriplovich and E. V. Pitjeva, Int. J. Mod. Phys. 
D 15, 615 (2006). 
[8] J. M. Frere, F. S. Ling and G. Vertongen, Phys. Rev. D 

77, 083005 (2008). 
[9] J. D. Anderson, P. A. Laing, E. L. Lau, A. S. Liu, 
M. M. Nieto and S. G. Turyshev, Phys. Rev. Lett. 81, 
2858 (1998). 

[10] S. L. Adler. larXiv:0805.2895l [astro-ph] . 

[11] E. Aprile et al, New Astron. Rev. 49, 289 (2005); 
Z. Ahmed et al. [CDMS Collaboration]. larXiv:0802.3530l 
[astro-ph]; R. Bernabei et al. [DAMA Collaboration], 
larXiv:0804.274TI [astro-ph], 

[12] F. J. Kerr and D. Lynden-Bell, Mon. Not. Roy. Astron. 
Soc. 221, 1023 (1986). 

[13] R. P. Oiling and M. R. Merrifield, Mon. Not. Roy. Astron. 
Soc. 297, 943 (1998). 

[14] J. F. Navarro, C. S. Frenk and S. D. M. White, Astro- 
phys. J. 462, 563 (1996); J. F. Navarro, C. S. Frenk and 
S. D. M. White, Astrophys. J. 490, 493 (1997); B. Moore, 
F. Governato, T. Quinn, J. Stadel and G. Lake, Astro- 
phys. J. 499, L5 (1998). 

[15] A. Klypin, H. Zhao and R. S. Somerville, Astrophys. J. 
573, 597 (2002). 

[16] M. S. Alenazi and P. Gondolo, Phys. Rev. D 74, 083518 
(2006). 

[17] M. D. Kislik, Kosmicheskie Issledovaniya 6, 853 (1964); 

V. V. Radzievskii, Sov. Astron. 11, 128 (1967). 
[18] M. Kamionkowski and S. M. Koushiappas, Phys. Rev. D 

77, 103509 (2008). 
[19] I. B. Khriplovich, Int. J. Mod. Phys. D 16, 1475 (2007). 
[20] J. I. Read, G. L ake, O. Agertz and V. P. Debattista, 

larXiv:0803.27T4l [astro-ph]; T. Bruch, J. Read, L. Baudis 

and G. Lake. larXiv:0804.2896l [astro-ph] . 



